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Abstract 



In this article the weU known "Perron-Frobenius theory" is investi- 
gated involving the higher rank numerical range ( A) of an irreducible 
and entrywise nonnegative matrix A and extending the notion of ele- 
. ments of maximum modulus in Ak{A). Further, an application of this 

' theory to the Afe(L(A)) of a Perron polynomial L{X) is elaborated via 

^ . its companion matrix Cl- 

> I 
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> 

oo ■ 1 Introduction 

(N 

^ ■ Let A^„(C) be the algebra of matrices A = [aij]^j^-^ with entries aij G C 

^ . and > 1 be a positive integer. The k-rank numerical range A/c(A) of a 

Q I matrix ^ G is defined by 

(1.1) Ak{A) = {A G C : PAP = XP for some PeVk}, 

where Vk is the set of all orthogonal projections P of C" onto any k- 
^ ', dimensional subspace /C of C". Equivalently, 

(1.2) Ak{A) = {X£C:X*AX = XIk,XeMnAQ^X*X = Ik}. 

For any k, the sets Ak{A) are generally called higher rank numerical range. 
The concept of higher rank numerical range has been introduced by Choi et 
aHn [21 [3l m [5] and studied thoroughly by other researchers in [111 IHl [IZ| • 
Apparently, for k = I, Ak{A) yields the classical numerical range of a matrix 
A 0|9], i.e. 

(1.3) Ai{A) = F{A) = {x*Ax : x G C", x*x = 1} 
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and it is readily verified 



F{A)^A2{A) 2 ... 5 Afe(^). 
Moreover, the notion of the numerical radius for the numerical range, 

r{A) = max{|2;| : z £ F{A)}, 

has been extended to the rank k-numerical radius 

rk{A) = max{|2;| : z £ Ak{A)}, 

considering rjt(^) = — oo, whenever Ai^{A) is an empty set [TT]. Note that 
r(A) > rk{A) and r{A) > p{A), where p(-) denotes the spectral radius of a 
matrix, i.e. p{-) = max{|A| : A € cr{-)}, with cr(-) to be the spectrum of a 
matrix. 

We mention that an n x n matrix A is said to be nonnegative when each 
aij > 0, and this is denoted by writing ^ > 0. Similarly, A is said to be 
positive whenever each aij > 0, denoted by A > 0. The matrix A G A4n is 
called reducible when there is a permutation matrix P such that 



P'^AP 



R S' 
T 



where R,T are both square. For n = 1, should be ^ = 0. Otherwise, A 
is said to be irreducible. If A is nonnegative and irreducible having q > 1 
eigenvalues of maximum modulus, then it is called imprimitive and q is 
referred to as index of imprimitivity. In case q = 1, yl is characterized as 
primitive. 

It is well known that Perron- Frobenius theory concerns the spectral prop- 
erties of positive and nonnegative matrices, namely the existence of positive 
or nonnegative eigenvalues and eigenvectors [8]. In addition, Issos' treat- 
ment contributes extensions of the Perron-Frobenius theorem to the numer- 
ical range of a nonnegative and irreducible matrix A, relating the p{A) with 
the r(A) [lOj. These results give the motivation for further investigation of 
Ak{A) in the case of a nonnegative and irreducible matrix A (section 2). 

In the next section, we present applications of the Perron-Probenius the- 
ory derived for a matrix polynomial, considering the higher rank numerical 
range of matrix polynomials. 

2 Nonnegative and irreducible matrices 

According to Issos' main theorems [TUl Th.4,7], if the numerical range F(A) 
of a nonnegative and irreducible matrix A has q maximal elements, then they 
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are equally spaced around a circle centered at the origin through a constant 
angle, with one of them lying on the positive real axis. In particular, 

j:{A) = |r(A)e'^ : t = 0, . . . , g - 1} 

is exactly the set of all the maximal elements in F{A) and to the numerical 
radius r{A) G F{A) there always corresponds a positive unit vector. Fur- 
thermore, Issos also proved that the cardinality of the set J-{A) coincides 
with the index of imprimitivity of A. 

In this section, we investigate to what extend do Issos' results apply to 
Afc(^) of a nonnegative and irreducible matrix A, taking into consideration 
that it is a non empty set At this point, we define the maximal elements 
in Afc(A) to constitute the set: 

Fk{A) = {ze Kk{A) : \z\ = rk{A)] , 

which for A; = 1 is equal to J- {A). Although the numerical range F[A) of 
a nonnegative matrix A always contains the numerical radius r[A) > (lUt 
Th.l], it is not generally true that rk{A) G Afc(A), as we may observe in the 
following example. 



"0030010 0' 
50000010 
00019006 

Example 2.1. Let the 8x8 nonnegative irreducible matrix A= 01000000 

00041002 
00100300 
.0500000 0. 

The outer curve in the following figure illustrates the boundary of F[A), 
whereas the second and third inner curves illustrate the boundary of A2(^) 
and A3(j4), respectively. Apparently, < r{A) € F{A), but in the figure, 
we easily recognize that < r2{A) ^ h.2{A) and < r:i{A) ^ ^.■i{A). Note 
that A has 4 maximal eigenvalues, which are marked by "+" and F[A) has 
4 maximal elements, as well. 




The following Lemma generalizes a familiar condition for rotational in- 
variance and symmetry about the origin to the case of the higher rank nu- 
merical range. 



3 



Lemma 2.2. Let A € A^„(C) be permutation similar (hence unitarily sim- 
ilar) to the matrix 



(2.1) 



C 



■ Ci2 
C23 



Cgl 










with the zero blocks along the main diagonal be square. For j = 1, . . . ,k such 
that rj{A) ^ {— oo,0}, we have: 

I. AjiA) = Ajie'^A), for t = 0,1,... ,q-l 

II. fi G Aj{A) if and only if fie 1 G Kj{A), for t = 0, 1, . . . , g — 1 

III. Aj(j4) is symmetric with respect to the origin, if q = 2t. 

Proof. I. By the proof of Theorem 6 in [lOj, D^^CD = e^^C, where D = 
/nj e e'^/„2 e . . . e e'(^"^)^I„^ is unitary diagonal with 9 = and ni + . . . + 



'"2 



Ua 



n, where ni {I = 1, ... ,q) is the dimension of the Z-th diagonal block of 
C . Hence, e^^ A = [P^ D P)~^ A{P^ D P) , where P is the permutation matrix 
such that A = P^CP, and Aj{A) = Aj{e''^A), for any j = 1, ... , k. 

II. By statement (I), we have that 

fi € Aj{A) ^ e'^/u G Aj{e'^A) = Aj{A), 
for alH = 0, 1, . . . , q — 1 and j = 1, . . . , k. 

III. Let q = 2t, then by (II), fi G Aj{A) if and only if e'^// G Aj{A), i.e. 



= -/i G Aj(^) for j = 1, . . . , /c. 
Example 2.3. Let the matrix A = 



□ 



as in ()2.ip with g = 3. The 



■ j 0" 

-j 
0002 
.3+2i 1 o_ 

boundary of A2{A) is illustrated below by the arched triangle. Apparently, 
A2{A) is rotationally invariant about the origin through an angle of ^ for 
t = 0, 1, 2, confirming Lemma l2.2l fI.II). 
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Any matrix A being permutation similar to the matrix C in ()2.ip is said 
to be q-cyclic and the largest positive integer q such that A is (7-cyclic is 
called the cyclic index of A. 

Regarding the number of maximal elements in {A) and their location 
on the complex plane, we refer to the following result. 

Proposition 2.4. Let A G A^„(]R) be imprimitive with index of imprimi- 
tivity q > 1 such that rk{A) > 0. Then 



(2.2) 

for every j 



2wt\ 

: t 



0, 



1 



1, 



, k with 9i = or 



Proof. Since 5 > 1 is the index of imprimitivity of A, there is a permutation 
matrix P such that 



P^AP 








Cql 



Cl2 

C23 








a 







Let 6j € [0, 27r) be the principal argument such that < rj{A) € Aj(e A) 



with j = 1, . . . ,k. Then, by Lemma [2.2r il). we have rj{A)e 
for t = 0, 1, . . . , g — 1 and j = 1, . . . ,k, whereupon we obtain 



ZTTt \ 



gA,{A) 



J^j{A)^\rjiA)e' 



t = 0,...,q-l 



for every j = 1, . . . , k. 

The index of imprimitivity q is equal to the largest positive integer such that 

;27r 

A is unitarily diagonally similar to the matrix e 1 A, equivalently, the ma- 



trices e ^^^A and e 



A are unitarily diagonally similar for the largest 



positive integer q. Therefore, the set 

2n 27r(g - 1) 
u, ) • ■ • ) 

q q 

is the cyclic group modulo 2it of the largest order, concluding that there 



does not exist = ^ < ^ such that rj{A)e^'i' G Kj{e~^^^ A) for j 
Hence we establish the equality 



1, 



^k. 



(2.3) 



J-^.(^) = |rj (A)e'(^^ + ^) : t = 0, . . . , g - 1} . 



In addition, if we denote by 
Kj{A) = Aj(yl), since A G Mn 



the conjugate of a set, it is clear that 
, i.e. Aj(A) is symmetric with respect to 



the real axis. Due to this symmetry and the equation (j2.3p . if we consider 
Oj 7^ 0, we obtain 27r — 9j = 6j + H!lli_Ji_ Hence, 9j = ^ and the proof is 
complete. □ 
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In view of the preceding proposition, the number of elements of maximum 
modulus in each Aj{A) of a nonnegative and irreducible matrix A, for j = 
1,2, ... ,k, is equal to the index of imprimitivity q. Especially, they are all 
successively distributed around a circle centered at the origin through the 
constant angle of Furthermore, it is also observed that after a clockwise 
rotation of Aj{A) about the origin through the angle of ^, it is achieved 

< rj{A) G Ajier'^A). 

On the other hand, if we consider the primitive class of nonnegative ma- 
trices A {q = 1), then the sets Aj{A), for j = 2, . . . , do not necessarily 
have only one maximal element, contrary to the F(A). This implication is 
justified by the symmetry of the sets Aj{A) with respect to the real axis, 
when rj{A)e^^i G Aj{A) for some 6j G (0, 27r), then also rj{A)e~^^^ G Aj{A). 



Example 2.5. Let the 8x8 imprimitive matrix A 



0020060 0^ 
10000070 
00023004 
3 000000 „rltVi 
3 000000 ^^"^ 
00040002 
10 3 
0900000 0. 

4. The boundaries of F{A) = Ai{A) D A2{A) D 



index of imprimitivity q 
A3(^) are illustrated below on the left figure and we verify that all Aj{A) 
have g = 4 maximal elements (j = 1, 2, 3), which are equally spaced through 
the angle of if, t = 0,l, 2, 3. Specifically, r2(A)e'5 G A2(^), rsiA) G A3(yl). 

^0 1 0" 
10 
10 
10 
1 
.1 1 1 0. 

with only one maximal eigenvalue. The boundaries of F{B) D A2(B) are 
illustrated on the right figure, revealing that A2 (B) has two symmetric maxi- 
mal elements, with respect to the negative real semi axis, whereas F{B) has 
only one maximal element r{B). The maximal eigenvalues of A and B are 
marked by "+" in both figures. 



For the primitive case, let the nonnegative irreducible matrix B = 





Remark 2.6. It is quite interesting to note that the n-cyclic permutation 
matrix P„ = ^"rT^l ^ -Aini^) is a special form of imprimitive matrix of 

t = 0, . . . ,n- 1. 
1}), where co(-) 



[0 Jn_l 

Li . 

index n, with spectrum the n-th roots of unity Zf = e^'^*'/", 
Apparently, P„ is unitary and F{Pn) = co({zo, z\,..., Zn- 
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denotes the convex hull of a set. By corollary 2.8 in 0, whenever 2k < n, 
we have that 

^k{Pn) = Co{{zo, Zi,..., Zn-l}), 

with zt {t = 0, . . . , n — 1) to be intersection points of the line segments 
[zt, zt+k] and [zt+i, zt+n-k+i], when Zj = zj-n for j > n - I. 
In addition, for the powers of Pn, we have that Ajt(P") = Afc(P„) (1 < a < 
n — 1), since the matrices are permutation similar to P„. Also, = In, 
whereupon Afc(P") = {!}. 

ro 1 0' 

For instance, F{P^) and A2(P5) for the 5-cyclic permutation P5 = 00010 

L 1 0. 

are illustrated below. 




Lemma 2.7. Let A G A1„(R) be nonnegative with irreducible hermitian 
part H{A) such that < rj{A) € Aj{e~'^^^A) for some dj G [0, 2tt) and every 
j = 1,2, . . . , k. If there exists an angle G M such that r{A)e'^ G P{^)j 
then rj{A)e^^ G kj{e~'^^ A) for every j = 1, . . . ,k. 

Proof. Suppose r{A)e^'^ G F{A) for some i;^ G M, then by proposition 3.7 
in [13] we have that e~^^A = D~^AD for some unitary diagonal matrix D. 
Obviously, < rj{A) G Kj{e-'^^ A) = Aj(e-'(^^+'^U), therefore rj(A)e^<^ G 
Aj{e~^^^ A) for every j = 1, . . . , A;. □ 

Investigating A^ [A) of a nonnegative matrix A in terms of the irreducibi- 
lity of the hermitian part H{A) of the matrix A, we extend an analogous 
discussion developed in [H]. 

Proposition 2.8. Let A G A > with H{A) be irreducible and 

rk{A) > 0. Then either J'j{A) coincides with the circle S{Q,rj{A)) for 

rj{A)e ^ 1 : t = 0, . . . ,q — 1> for every 

j = 1, . . . ,k, where 9j = or 6j = ^ and q > 1 is the largest positive integer 

i— 

such that A is diagonally similar to the matrix e 1 A. 

Proof. Suppose that for each j = 1, . . . , A; we have < rj{A) G Aj{e~^^^ A) 
for some principal argument 6j G [0,27r). Due to proposition 3.11 in [13] 
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; 27rt 
<3 



and corollary 3.6 in [14], either F{A) = 5(0, r(A)) or T{A) = |r(yl)e'' 

t = 0, . . . , g — 1}, where q is the maximum positive integer such that A is 

unitarily diagonally similar to the matrix e i A. The first equality of the 
sets indicates that r(A)e"^ G ^'{A) for every angle </) G M, thus by Lemma 
[2T71 rj{A)e^'t' G Aj(e~'^-'.4) for every angle (/> G M, concluding J^j{A) = 

5(0, rj(74)) for all j = 1, . . . , k. For the second case where r{A)e « G -^(^4), 

Lemma l2 . 71 verifies rj{A)e i G Aj(e~^ ^ A) and then 

{r,(A)e^('^+f ^ t = 0, . . . , g - l} C J-,(A). 

The equality of the sets, with 6j = or ^ is established similarly as in the 
proof of Proposition 12.41 since q is identified with the largest positive integer 

J/ Q . I 27r \ 

such that e i A is unitarily diagonally similar to the matrix e ^ « M. and 
kj{A) are symmetric with respect to the real axis for every j = 1, . . . ,k. □ 

Corollary 2.9. Let A G A^n(I^)) ^ > with irreducible hermitian part. If 
< rj{A) ^ kj{A) for some j = 2, . . . ,k and A.j{A) is not a circular disc, 

then kj{A) is symmetric with respect to the lines C± = |ze^*9 : z G m|, 
where q is the largest positive integer such that A is diagonally similar to 
e 1 A. 

( .7r(2t + l) ^ 

Proof. By Proposition l2.81 Tj{A) = < rj{A)e i : t = 0, . . . ,q — 1> , where- 
upon 

J-j(e^'f ^) = |rj(^)e'^ : t = 0, . . . , g - l} C 9Aj (e^'f A), 
for some j = 2, . . . , A:. Clearly, Aj(e ^ A) is symmetric with respect to M, 

qzj— dii — 

implying Aj{A) = e i Aj{e i A) to be symmetric with respect to the lines 
C±. □ 

Corollary 2.10. Let A G 7W„(M), ^4 > with H{A) be irreducible such that 
rk{A) > 0. If F{A) is a circular disc, then Aj{A) is also a circular disc for 
every j = 2, . . . ,k. 

Any matrix of the form ()2.ip with Cqi = and g > 1 is called to be 
a block-shift matrix. It is already known that for real nonnegative matrices 
with irreducible hermitian part, which can be put into the block-shift form by 
means of a permutation, we obtain the circularity of the numerical range 
Th.l]. Because of Corollarv 12.101 it is immediate that this type of matrices 
also characterize the circularity of the higher rank numerical range. 

Proposition 2.11. Let A G A1„(M), A>0 with H{A) be irreducible. If A 
is permutation similar to a block-shift matrix, then Aj(A) is identified with 
the circular disc T>{0, rj{A)) for j = 1, . . . , k. 
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We should note that the conclusion of the preceding proposition does 
not hold in general for any complex matrix A being unitarily similar to a 
block-shift matrix. For the general case, we present the next result. 



Ai 




with 



Proposition 2.12. Let A G A^„(C) be unitarily similar to 
Ai G A4m,n-m(C) and rankAi = k. Then 

A,(A)=P(0,^^)/ori = l,...,fc, 

where aj{Ai) denotes the j-th largest singular value of Ai and then rj{A) 
^,forj = l,...,k. 



Proof. It only suffices to prove that Aj{ 



Ai 




p(0,^), for j 



1, . . . ,k, since Aj(A) is invariant under unitary equivalence. Therefore, by 



Ai 




)= ^ e-'" ■{zeC:Rez<Xj{H{ 

6»G[0,27r) 



e'Mi 




where Xj{H{A)) denotes the j-th largest eigenvalue of the hermitian part 
H{A) of matrix A. Hence 



A,( 



Ai 




ee[o,27r) 



z G C : Rez < -Xj{ 



e'Mi 
e-^^Al 



It is known that the eigenvalues of the hermitian matrix 







e-^^Al 



are the singular values cJi(yli) > . . . > ak{Ai) > > . . . > > —ak{Ai) > 
. . . > —ai{Ai) with k = rankAi [8], thus for j = 1, . . . ,k 



A, 



Ai 




) 



n 



e-'^<!zGC:Rez<^^ 



ee[o,27r) 



P(0,^). 



□ 



3 Application to matrix polynomials 

An extension of the aforementioned results to the higher rank numerical 
range of a matrix polynomial arises naturally and it is the purpose of this 
section. For this reason, we refer to Perron polynomials L{X), which are 
n X n monic matrix polynomials of mth degree 

(3.1) L(A) = /A'" - Am-iX"'-' -...-AiX-Ao, 
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with Aj to be nonnegative matrices, for j = 0, . . . ,m — 1. The higher rank 
numerical range Afc(L(A)) of L(X) has been recently defined in [1] by the set 



(3.2) Afc(L(A)) = {A G C : Q*L{X)Q = 0^ for some Q G Mn,k, Q*Q = 4} , 
which for k = 1 yields the numerical range 

(3.3) w{L{X)) = {A G C : x*L{X)x = for some x G C", x*x = 1} . 

The notion of the Perron polynomial L{X) in ()3.ip is equivalent to the non- 
negativity of its mn x mn companion matrix 



Cl 



/„ 
/„ 



^0 








A, 



m—l 



hence an extension of the main Perron-Frobenius theorem concerning the 
spectrum a{L) of L{X) and Issos' results on the numerical range ■w{L{X)) 
are established via the companion matrix Cl jl5j . Further, we denote the 
rank /c-numerical radius of Afc(L(A)) 

rfc(L) = max{|A| : AG Afc(L(A))}, 

with rk{L) = —oo when Afc(L(A)) = and consequently, the set of maximal 
elements in Afc(L(A)): 

7-,.(L) = {A G Afc(L(A)) : \X\=rk{L)}. 

As we have noticed in section 2, clearly Afc(L(A)) does not always contain 
the element rk{L) > 0. 

The following lemma is a generalization of a corresponding result in [15] , 
which identifies Afc(L(A)) with a specific subset of Afc(Ci) of the companion 
matrix Cl of a matrix polynomial L{X). We should recall that Ak{CL) 
always contains Afc(L(A)) as it has been proved in [H Prop. 16]. 

Lemma 3.1. Let the mn x k matrix Y{X, Q) = (1, A, . . . , X"''^)'^ (g) Q with 
A G C and Q G Mn,k- Then 



Afc(L(A))U{0} = {X:Y*{X, Q)ClY{X, Q) = A4, Q G Mn,k , Q*Q 



c(|A 



27 -^fc 



where c{\X\^) = (1 + \X\^ + \Xf + ... + |Ap("-^V/2. 
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Proof. It is readily verified that 

y*(A, Q){Cl - Mmn)Y{X, Q) = ~X'^-'Q*L{X)Q 

a nd Y*{X ,Q)Y{X,Q) = c{\X\'^ )Q*Q. Hence, Aq G Afc(L(A)) if and only if 

\J c{\Xo\'^)Q* L{Xo)Q\J c(|Ao|^) = 0^ and equivalently when Q*L{Xq)Q = Ofc. 

□ 

Proposition 3.2. Let L{X) be a Perron polynomial as in (j3.ip such that 
rk{L) > 0, with irreducible companion matrix Cl and q > I eigenvalues of 
maximum modulus. Then the maximal elements in Aj{L{X)) are exactly the 
set 

(3.4) T,{L) = [r,{L)e<'^^"-f\ t = 0, . . . ,q - l] , 
where 6j = or 9j = ^, for j = 1, . . . ,k. 

Proof. The companion matrix C^, > is an imprimitive matrix with in- 
dex of imprimitivity q > 1, since cr(L) = a{CL). Hence q is the largest 

;27r 

positive integer such that Cl is unitarily diagonal to e i Cl- Therefore, 

i— _i 

e 1 Cl = D C^D (t = 0, 1, . . . , g' — 1), for some unitary diagonal mn x mn 

. 27rt ■ 27rt(m-l) 

matrix D = diag{D(), e i Dq, . . . , e i Dq), with Dq G A4n to be uni- 
tary diagonal. Suppose Oj S [0, 27r) is the principal argument such that 
< rj{L) G Aj(e"'^^L(A)), then by Lemma [3TT] we have 

Ok = Y*{rj{L)e''^,Q){CL - rj{L)e''^ I^n)Y{r,{L)e''^ ,Q) 

= Y*{rj{L)e''^,Q)D-\e-''-fCL - r,{L)e''^ Imn)DYirj{L)e''\Q) 

= Y*{r,{L)e'^'^^"-f\DoQ){CL - rj{L)e'^'^^fhmn)Y{r,{L)e'^'^^"-f\DoQ), 

with Y*{e'''fx,DoQ)Y{e''fx,DoQ) = Y*iX,Q)Y{X,Q) = Imn,k for any 
scalar A. Apparently, 

(3.5) Fj{L) D [rj{L)S''^'-f\t = 0,...,q-l]. 

Further, as in the proof of Proposition 12.41 we establish the equality of the 
sets with Oj = or 0j = ^ for j = 1, . . . , /c, since the coefficients of L{X) are 
real, whereupon Afc(L(A)) is symmetric with respect to the real axis. □ 
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